Strongly Intensive Quantities 



M. I. Gorenstein^'^ and M. Gazdzicki^' ^ 

^ Bogolyuhov Institute for Theoretical Physics, Kiev, Ukraine 
^Frankfurt Institute for Advanced Studies, Frankfurt, Germany 
^Institut fiir Kernphysik, University of Frankfurt, Frankfurt, Germany 
"^Jan Kochanowski University, Kielce, Poland 

Abstract 

Analysis of fluctuations of hadron production properties in collisions of relativistic particles profits 
from use of measurable intensive quantities which are independent of system size variations. The first 
family of such quantities was proposed already in 1992. The second is introduced in this paper. We 
also present a proof of independence of volume fluctuations for quantities from both families within 
the framework of the grand canonical ensemble. These quantities are referred to as strongly intensive 
ones. Influence of conservation laws and resonance decays is also discussed. 
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I. INTRODUCTION 



An intensive quantity is a physical quantity which does not depend on the system volume. 
By contrast, an extensive quantity is proportional to the system volume. Clearly, the ratio of 
two extensive quantities is an intensive one. For example the number of particles, A^, in the 
relativistic gas fluctuates around its mean value, (A^). Within the grand canonical ensemble (A^) 
is an extensive quantity, whereas the ratio of mean multiplicities of two different particle types 
is an intensive one. Particle number fluctuations are quantified by the variance, (A^^) — (A^)^, 
which is an extensive quantity, while the scaled variance, [(A^^) — (A^)^]/(A^), is an intensive 
one. 

Statistical models are surprisingly successful in modeling multi-particle production in high 
energy interactions They are used to describe properties of strongly interacting matter 
created in nucleus-nucleus collisions in terms of intensive quantities. In particular, an equation 
of state is usually given as a function relating pressure, temperature and baryonic chemical 
potential. On the other hand, in high energy collisions the volume of produced matter cannot 
be kept fixed. For instance, nucleus-nucleus collisions with different centralities may produce 
a statistical system with the same local properties (e.g., the same temperature and baryonic 
chemical potential) but with the system volume changing significantly from interaction to inter- 
action. Thus, an important question is whether it is possible to measure the properties of the 
system without knowing its volume fluctuations, or equivalently, whether there are measurable 
quantities which are independent of volume fluctuations. 

Within the grand canonical ensemble the answer is yes, and the quantities with the required 
properties will be referred to as strongly intensive ones. Ratios of mean particle multiplicities 
are strongly intensive quantities. In general, this is the case for all ratios of any two extensive 
quantities which correspond to the first moments of fluctuating variables. They are intensive 
and strongly intensive quantities. The situation is, however, more complicated for the measures 
of fluctuations which include the second moments of fluctuating variables. For example, as will 
be shown below, the scaled variance of a particle number distribution is an intensive quantity, 
but not a strongly intensive one. 

In this paper we show that there are two families of strongly intensive quantities which 
characterize the second moments of random extensive variables used to study fluctuations and 

n 

correlations in a physical system. While the first family was introduced already in 1992 [2[, the 
second one is proposed in this paper. 

The paper is organized as follows. In Sec. II the two families of strongly intensive quantities 
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are introduced. For simplicity, this is done within the model of independent particle sources. 
The relation of the strongly intensive quatities to previously used fluctuation measures is dis- 
cussed in Sec. III. The proof that the quantities are in fact strongly intensive, i.e., strictly 
independent of volume and volume fluctuations within the grand canonical ensemble is given 
in Sec. IV. Finally, their properties within the canonical and micro-canonical ensembles are 
discussed in Sec. V. Summary given in Sec. VI closes the paper. 



II. TWO FAMILIES OF STRONGLY INTENSIVE QUANTITIES 

Let us start from the model of independent sources for multi-particle production in which 
the number of sources, Ng, changes from event to event. In this model, extensive quantities 
(e.g., mean number of particles, mean transverse energy) will be considered as those which are 
proportional to A'^s- Two fluctuating extensive variables A and B can be expressed as: 

A = ai + a2 + . . . + aN, , B = bi + h2 + . . . + h^, ■, (1) 

where and hk denote the contributions from the /c-th source. One finds for event averages: 

{A) = (a) (AT.) , {A') = {a') (iV.) + {af [(iVf) - (iV.)] , (2) 
{B) = {b) (iV.) , {B') = {b') (iV.) + {bf [{N^} - (iV.)] , (3) 
{AB) = {ab) (iV.) + (a)(6) [(iVf) - (iV.)] , (4) 

where (a), (6) and (a^), (6^), (ab) are the first and second moments of the distribution P*{a,b) 
for a single source. These quantities are independent of A^^ and play the role of intensive 
quantities in the model of independent sources. The distribution P*(a, b) is assumed to be the 
same for all sources, i.e., they are statistically identical. The probability distribution Ps{Ns) of 
the source number is needed to calculate (A''^) and (A^^ ) and, in general, it is unknown. Using 
Eq. ([2]) the scaled variance ua which describes the event-by-event fluctuations of the extensive 
variable A can be presented as: 

u^A = = + (a) = + (a) c.. , (5) 

where u* is the scaled variance of the quantity A for each source. A similar expression follows 
from Eq. for the scaled variance ub- The scaled variances ua and ub are independent of 
the average number of sources (Ng). Thus, oja and ub are intensive quantities. However, they 
depend on the fluctuations of the number of sources via Us and, therefore, they are not strongly 
intensive quantities. 



From Eq. (|2]) follows that a knowledge of (A) and (A^) is not sufficient to derive any strongly 
intensive quantity. This is, however, possible when two extensive random variables, A and B are 
considered. In order to characterize fluctuations, one may then construct special combinations 
of the second moments in which the terms proportional to (N^) in the r.h.s. of Eqs. 
are not present. Clearly, only two linearly independent combinations of this type result from 
three equations (EJH]). Note that in order to remove the dependence on (Ns), strongly intensive 
quantities should be in a form of reducible fractions. The following combinations seem the most 
convenient: 

= {C)-' [ (B) UA + {A)coB - 2i {AB) - {A){B) ) ] , (6) 
A^^ = (C)-' [ (B) UA - {A) UB ] , (7) 

where (C) is the average of any extensive quantity, e.g., {A) or {B). Straightforward calculation 
of dH]) and dlD using Eqs. (015]) gives: 

S^^ = (c)-^ [ (6) ujl + (a) ujl - 2{ {ah) - {a){b) ) ] , (8) 
A^^ = {c)-' [ {b) co: - {a) ujl ] . (9) 

Thus, S^^ and A^^ defined by Eqs. ([MZD depend on cj* and w^, but they are independent of 
the average number of sources {Ns) and its fiuctuations ujs- They are, in fact, strongly intensive 
measures which quantify fiuctuations of any two extensive random variables A and B. This will 
be proved in Sec. IV within the grand canonical ensemble for a case of particle multiplicities. 

There is an important difference between the S"^^ and A"^^ quantities. Namely, in order 
to calculate A"^^ one needs to measure only the first two moments: (A), {B) and {A^), {B"^). 
This can be done by independent measurements of the distributions Pa{A) and Pb{B). Quan- 
tity S"^^ includes the correlations term, {AB) — {A){B), and thus it requires, in addition, 
simultaneous measurements of A and B in order to obtain the joint probability distribution 
Pab{A, B). The quantities S"^^ and A"^^ also have different symmetry properties: T^^^ = T,^^ 
and A"^-^ = — A-^"^. We call all strongly intensive quantities which include the correlation term 
the E family, those which include only variances the A family. 

III. RELATION TO OTHER FLUCTUATION MEASURES 

The well-known fluctuation measure $ was introduced in 1992 [2;] for a study of transverse 
momentum fluctuations and it belongs to the S family. In the general case (sl, when A = X = 
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Xlili represents any motional extensive variable as a sum of single particle variables, and 



i? = is the particle multiplicity, one gets: 



(AT) 



1/2 



— 



1/2 



(10) 



where is given by Eq. ([6]) with C = N ^ and x^, correspond to single-particle inclusive 
averages. Note that these inclusive quantities can also be presented in terms of event averages, 
namely x = (.X)/{N) and = (X2)/(iV), where X2 = T^ti^l The measure $ was extended 
in 1999 {4, 5| for multiplicity fluctuations. For two particle types, A and B, the $ measure 
was constructed by setting Xj = 1 if the i-th particle is of the A type and Xi = otherwise. 
One then finds x = x^ = (A) /[{A) + (B)] and thus li^ - x^ = {A){B)/[{A) + {B)]^- with A 
and B denoting particle numbers. Taking into account these relations, and using X = A and 
N = A + B in Eq. (fTOj) . the expression for the $ reads: 



(11) 



where T,^^ is given by Eq. ([6]) with C = A + B. 

A possible extension of $ for the case of two motional variables was not discussed up to now, 
however it can be naturally done within the framework of the T,^^ and A^^ families presented 
here. It is also important to note that the $ measure extended to the study of the third moment 
preserves its strongly intensive properties within the model of independent sources [g^. Study 
of strongly intensive quantities which include 3''°' and higher moments of extensive quantities 
is beyond the scope of this paper. 

Another quantity frequently used to characterize the fluctuations of particle numbers A and 
B was introduced in 2002 [3] as: 



.AB 



(AiA-l)) , (fi(fi-l)) 



Using Eq. 



dyn ^^^2 

with C = A + B, one easily finds the relation: 

(A + B) 



(AB) 
{A){B) 



;i2) 



1^: 



AB 

dyn 



[ T.^^ - 1] 



(13) 



{A){B) 

Equation (IT^ shows that , similar to S"^-^, is independent of fluctuations of the source 
number, but it decreases as z/^^ oc {Ns)~^ and, thus, it is not an intensive quantity. Note that 
the quantity {C)^^^, where C can be chosen as A, B, or A + B, is a strongly intensive quantity 
from the S family. Despite the fact that specific examples of the T,^^ family were introduced 
and discussed a long time ago, the A"^^ family is proposed in this paper for the first time. 



IV. PROOF WITHIN THE GRAND CANONICAL ENSEMBLE 



Let us now prove within the grand canonical ensemble (GCE) that the two families of 
quantities, T,^^ and A^^ , are strongly intensive. The proof will be limited to a case of particle 
multiplicities, i.e., A and B will be number of particles of type A and B, respectively. The GCE 
partition function H of the quantum gas which is a mixture of different types of particles reads: 



expj V J2^J ^^^^ ^3 ^M-^j/T)] 



(14) 



where V and T denote, respectively, the system volume and temperature, Xj is the fugacity 
which is related to particle chemical potential as A-,- = exp(/ij/T), dj denotes the number of 
a particle internal degrees of freedom, ej = (m^ + p'^Y^'^ is the particle energy with rrij and p 
being its mass and momentum, rji = —1 for bosons, rji = 1 for fermions, and r]i = corresponds 
to the classical Boltzmann approximation. The GCE averages are calculated as: 

a=iA.^5^W^ /; ^V,u, (15) 

- CAa J [27iy A^ exp(eA/T) + r]A 

^=4 (>^a4-Xs = V^n\ + V dA-/eMeA/T) ^^^^ 



^^a) ^ J (27r)3 [A:,^exp(6Vr) + r^^]^ ' 



d , d ^ 



AB = - A^-— As-— = V' UAUB , (17) 
^ OXa oXb 

where ua = A/V and ub = B/V denote the particle number densities. The corresponding 
expression for B and B'^ are obtained by replacing A by 5 in Eqs. fll5m6p . 
For the GCE scaled variance one finds, 

~^ - a' f d^p dAX^'exp{eA/T) 



uja = = = n 



A " J (27r)3[A-iexp(6VT) + VaY ' 

It corresponds to the particle number fluctuations at a fixed volume V . It is an intensive 
quantity and depends only on T and ^ia- Note that o;^ > 1 for bosons, uj\ < 1 for fermions, 
and u\ = 1 for classical Boltzmann particles. 

We introduce now volume fluctuations assuming that local properties of the system within 
the GCE, i.e., the temperature and chemical potentials, are volume independent. The volume 
fluctuations will be described by the probability density function F(y). Thus, the full averaging 
denoted as (. . .) will include both the GCE averaging f lT5lfT7|l at a fixed volume and an averaging 
over the volume fluctuations: 

(A) = (V) UA , (A^) = {V^) n\ + {V) HA uj*A , {AB) = {V'') haUb , (19) 
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where {V^) = J dV V''F{V) for k = 1,2. One finds, 



uja 



^A + Ua UJv ■ 



(20) 



(.1) ' - {V) 

The corresponding expression for is obtained by replacing A by 5 in Eq. (120|) . Equations (120|) 
and ([5]) have a similar structure. Namely, the first terms or correspond to the particle 
number fiuctuations at a fixed volume V or fixed number of sources Ns, respectively. The 
second terms correspond to the contribution from the volume fiuctuations in Eq. fl20|) and the 
fiuctuations of the number of sources in Eq. (E]). 

Calculating (EHZ]) according to Eq. f lTI?]) with C = A + B one gets: 



1 



ua + riB 



Ub u*a + riA UJ*B 



A 



AB 



Ub i^A - ^A ^B 



(21) 



Ha + Ub 

Equation fl2T]) proves that S^"^ and A"^^ are strongly intensive quantities as they are strictly 
independent of average volume (V) and its fiuctuations ujy- 

Note that in the GCE there are no correlations between the number of different particle 
species, i.e., {AB) — {A){B) = 0. In the applications to hadron production in high energy 
collisions stable particles are detected, whereas the GCE system includes also short lived reso- 
nances which finally decay into stable particles. These resonance decays increase multiplicities 
of stable particles and thus change numerical values of ua, tib, w^, and Ub- If decay products 
of a resonance R decay include both A and B hadrons a correlation between them appears and 
it can be expressed as: 



{AB) - {A){B) = Y,{R)\{AB)n- {A)n{B)n\ ^ J^^^) P% 



R 



(22) 



R 



where {R) is a mean multiplicity of R and {■ ■ ■)r means the averaging over the its all decay 
channels. The measure T,^^ will have then the form: 



^AB 



Hb UJa + ^. 



Ha + Ub 

and thus it remains a strongly intensive quantity. 



B 



nR Pab 



R 



(23) 



V. PROPERTIES OF S^^ AND A^^ WITHIN CANONICAL AND MICRO- 
CANONICAL ENSEMBLES 

For a large volume system in equilibrium the particle number distribution P{A, B; V) in 
the GCE, canonical ensemble (CE), and micro-canonical ensemble (MCE) can be written in a 
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general form of the bi-variate normal distribution (see Refs. 
P{A,B-V) = ^ [col co% {1 - p'X^) AB\-''' 



0,3): 



X exp 



1 



{A^Af {A-A){B-B) 
^ ^Pab 



[B-B) 
u%B 



(24) 



2 (1 - P*L) \ co*aA " K ^*B A Byi^ 
where A = haV and B = ubV are mean particle numbers. Averaging at a fixed volume V 
is defined as [k = 1, 2): 

A^ =J2^''PiAB;V) , B^ =J2b''P{A,B;V) , AB = ^ P(A, S; 1/) . (25) 



A,B 



A,B 



A,B 



The straightforward calculations of f l25p with the distribution function f l24p give: 



A 



52 



B 



AB - AB 



Pab 



A ' 'B ^ r * * ^"r5ii/9 ~ Pab ■ (26) 

A B Pa^s ^ -BJ / 

Equation (1261) reveals the meaning of the parameters in the distribution - the scaled vari- 
ances u\ and cuIj, and the correlation coefficient p^^. In a mixture of relativistic ideal gases 
particle numbers are not conserved, and thus A and B fluctuate in all statistical ensembles. 
This leads to non-zero positive values of oj^ and which are approaching constant values 
with system volume increasing to infinity^. Particle correlations, and thus non-zero p^^, result 



from exact material and motional conservation laws Thus the correlation coefficient 
equals to zero in the GCE, and is non-zero in the CE and MCE. The exact conservation laws 
also influence the values of and co^. The quantities like the particle number densities do 
not depend on the choice of the statistical ensemble for large systems. This means thermody- 
namical equivalence of the statistical ensembles. Below we present results for the CE and MCE 
in the large volume limit in which all three statistical ensembles become thermodynamically 
equivalent. Let us however stress that the thermodynamical limits of the quantities fl2B|) are 
different (see Ref. [lli for details) in the GCE, CE, and MCE ensembles. 

We introduce now the volume fluctuations assuming that local properties of the system (e.g., 
temperature and conserved charge densities in the CE, or energy density and conserved charge 
densities in the MCE) are volume independent for large enough volumes. In this case, the 
distribution f l2^ depends on the system volume only through the average multiplicities. The 
full averaging reads: 



(...) 



f dVFiV) PiAB;V) . 

A,B 



(27) 



^ The volume dependence may be different for system at the phase transition. For example, in a case of the 
Bose-Einstein condensation one gets, uj oc 1/1/3 at T = Tc and w oc F at T < Tc, as shown in Ref. 0. 
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Calculating ([MZl) according to Eq. fl27|) with C = A + B one gets: 

S^^ = 7 riB u;\ + - 2p\^ {uaUb w^^b)^^^ > (28) 

A^^ = 1 [nBUj\ - nAUj*B] ■ (29) 

Equations f l28p and (12^ show that S"^^ and A"^^ are independent of the average system volume 
and its fluctuations. For the CE and MCE this is valid if the volume fluctuates in a range in 
which all three statistical ensembles are thermodynamically equivalent. 

A unique determination of flve intensive quantities, tt-a, n^, oj% u^, and p\^, from measure- 
ments of {A), {B), {A"^), {B^), and {AB), is impossible as {V) and uy are in general unknown. 
The average particle multiplicities are given by (A) = nA{V) and (B) = nB{V). Therefore, 
only the ratio of particle number densities, tab = {^) / (B) = ua/tib, can be found from the 
measurements of {A) and {B). The three strongly intensive quantities, tab, S^"^, and A^^, 
allow a unique determination of u\ and in the GCE, if resonance decay effects are absent. 
In this case there are no correlation between A and B at flxed volume, and from Eq. fl2T|) one 
flnds: 



2 (1 + tab) [S^^ + A^^] , co*^ = 2{1 + t^'b) [S^^ - A^^] . (30) 



However, if p*^^ ^ 0, as in the CE and MCE, or including correlations due to resonance decays in 
the GCE, even the knowledge of all strongly intensive quantities is not sufficient to reconstruct 
u\, Ub, and Pas ^ unique way. 

VI. SUMMARY 

In summary, in this paper we consider two families of strongly intensive quantities T,^^ and 
A^"^ which characterize fluctuations of system properties. While speciflc measures from the 
S"^-^ family were introduced already in 1992 [2], the A^"^ family is proposed in this paper for the 
flrst time. We prove within the grand canonical ensemble that both T,^^ and A"^^ quantities are 
strictly independent of volume and volume fluctuations. In the canonical and micro-canonical 
ensembles they are approximately independent of volume and volume fluctuations for large 
enough systems. Note that the $ and measures, which can be expressed in terms of 
T,^^, have already been used successfully to study transverse momentum and particle ratio 



fluctuations (see, e.g., Refs. [l2|, ll3[ and references therein). We hope that the results presented 
in this paper will be useful in further analysis of fluctuations of hadron production properties 
in collisions of relativistic particles. 
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